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ABSTRACT 

[UnclaBsified] 

Thin, layered, and semi'inflnlte targets have been subjected to high velocity 
projectile impacts experimentally, and theoretical calculations performed by the 
DORP code for these same target impacts. The DORP code is a two-material, 
two-dimensional continuous Eulerlan hydrodynamic code coupled vrith a rigid, 
perfectly plastic strength model. In addition, DORFO, a nine-material version of 
DORP, has calculated several of the impact experiments. The agreement is good 
between the theoretical calculations and experimental values for shock attenu¬ 
ation and projectile length loss. 


PROBLEM STATUS 

This is a final report on this phase of the project; work is continuing on 
other phases. 


AUTHORIZATION 

NRL Problem F04-21 
ARPA Order 864, Program Code 8E20 
Project N00014-70-C-0002 


Manvucrlpt tubmiUed February 1, 1972. 




THEORETICAL AND EXPERIMENTAL IMPACT STUDIES 
[Unclassified Title] 


INTRODUCTION 

Significant improvements have been made in the strength and transport formulation 
of the DORF code (1). In addition, the code has been modified to allow the treatment 
of up to nine different equations of state in the same problem (2). 

It was felt that at some point experiments should be performed to substantiate the 
results from the code calculations. As a result, several impacts into thin, layered, and 
semi'infinite targets have been performed experimentally at NRL. Shock-propagation, 
shock-attenuation, and projectile-loss compa^ons between theory and experiment have 
shown very good agreement. 

As a result, we feel confident that numerical techniques, such as the DORF code, 
can be applied to a variety of interesting problems with a high degree of success. 


THE HYDRODYNAMICS AND STRENGTH 
FORMULATION FOR THE DORF CODE 

The DORF and DORF9 codes numerically solve the Eulerian equations of hydro- 
dsmamic flow, which are basically statements of mass, momentum, and energy conser¬ 
vation respectively: 

V • (pu)-O 

~ + V • (PMU)--V P 

OX 

boE 

+ V ' ipEit) •> — V • (Pu), 

where 

u ■ flux vector 
p ■> material density 
P * hydrostatic pressure 
E ■ specific internal energy. 

As a separate step, the momentum and energy conservation equations for material 
strength are solved: 
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3U,. 

P 17 “ OiH 
bE 

^ Ti" 

where Oj^ U the stresa deviation tensor and the summation and implied differentiation 
conventions apply. An elementary flow chart showing the manner and order in which 
DORF and DORF9 solve these equations is given in Appendix A. Reference 1 gives a 
more complete description of the methodology. 


THE EQUATION-OF-STATB MODEL 

USED IN THE DORF AND DORF9 CODES 

The equation of state used in the DORF and DORF9 codes is given in a report by 
Tillotson (3); it yields pressure F as a function of volume V and specific internal energy B. 
For condensed states, when p/pq > 1 or for any cold state when B < E,, the equation 
becomes 


P “ Pc ■ Ep no +-^-\ + Au + Bp®, 

\ E/B^n^ + IJ 

where 

Po * initial (zero pressure) density of material 

E, " energy required to bring material to vaporization point 

V “ P/Po 
p - rj — 1 

and a, b. A, B, and Eq are empirical fits to the experimental data. A cutoff is built into 
the codes which does not allow the use of p < — A/2B in this equation since this could 
allow 9P/3V to be positive. 

For expanded states, where p/pq < 1 and £ > £,’, the equation of state is 

, r ^^P . . ~ 111 ~ 11* 

P'-Ps “aEp + ^--- + Ape e " , 

l(E/Eon2) + 1 J 

where B^' - E^ plus the energy of vaporization, and a and 0 are empirical constants. 

In the intermediate region, where r; < 1 and E^ < E < E, & smooth transition 
between the condensed and expanded states is assured by setl< g 
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When the term E/Eqj]^ « 1, the equation for Pq behaves like a Mie-Gruneisen 
equation of state vnth a constant Qruneisen ratio of a + b. At large energies and mod¬ 
erate compressions, E/EqT)^ » b and n small, Pc behaves like a gaseous equation of state 

P “ apE, 

The Gruneisen ratio over ali energy and densities for which P is valid is given by 

Q{E, p)s vl^] ~a + - - -. 

\a^/v («/JEo»?2)+l 

The yieid strength is represented by 

where a' is approximately 0.07 for metals and E^ is the energy required to melt the 
material, If E is greater than E„, V is set to zero. 

The stress deviators are defined as 

Oij - Biij, 

where 

T/iipTj 

Ujj • stress deviation tensor 
ijj “ strain rate deviation tensor 

and the summation convention again applies. 


COMPARISON BETWEEN CALCULAHONS AND EXPERIMENTS 

Shock Propagation into Semi-Inflnite Targets for Like and 
Unlike Material Hypervelocity Impacts (DORF and DORP9) 

Two shock propagation experiments, hereafter referred to as shot 1 and shot 2, 
were performed, and the corresponding code calculations completed using DORF and 
DORF9. 

Shot 1 


A .3/8-in.-diameter aluminum sphere was impacted normally into a semi-infinite 
HOOF aluminum target at 7.091 km/sec. The decay of the on-axis peak pressure of 
the shock wave was observed using manganin pressure gauges imbedded at various points 
within the aluminum target (4). Additional experimental data was obtained from the 
work of Charest (6), who obtained peak pressures using free-surface particle velocity 
techniques. Both sets of data are shown in Fig. 1. The DORF code was then used to 
calculate this same impact, and t' n calculated nressures were compared with the 
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experimental values. The results were somewhat disappointing at first because the DORP 
code values were consistently lower than either of the experimental measurements, although 
there was a convergence In the lower pressure regions (large distances from the impact 
point). Reasons for the discrepancy were difficult to find, and suspicion finally settled 
on the Tillotson equation of state for aluminum. 

A plot of shock velocity t/, versus particle velocity Up was derived from the Tillotson 
formulation, and the results were compared witli the experimental data of Rice, McQueen, 
and Walsh (6) and of Al’tshuler et al (7). As shown in Fig, 2, the Tillotson equation 
starts deviating from the experimental values fairly quickly, reaching a value for U, that 
is about 8% low at the Hugoniot conditions for shot 1 {Up »» 3.6 km/sec.). This means 
that the Tillotson equation of state would predict a Hugoniot pressure that is 8% low 
for shot 1. This, combined with the generally more compressible (than reality) nature 
of the Tillotson formulation over the whole range of compression, could conceivably 
affect both the magnitude and decay rate of the calculated shock wave. To determine 
whether or not this reasoning applied, a simple Mie^Gruneisen equation-of-state option, 
using the Hugoniot as the reference pressure curve, was Inserted into DORF. The 
equation of state used the following standard equations: 
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+ Gp(E-Eh) 


p . ^0^0^ 

^ ’ (l-AW)^' 
■e 

Efi - (l/2)f>«(Vo - V) 

“ 1 - ifio/P) 

Po ■ 2.71 g/cc 
O - 2.13 
Cq ■> 6.38 km/sec 
A • 1.337. 
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Fig, 3 — OORF code calculaMoni vi 
eKperimenUl data for shot ] 


The laet two conatants come from 


(/, - Co + AUp 

and fit the experimental data quite well as can be seen in Fig. li. 

Since our Impact is basically a shock propagation problem, the material values near 
the region of interest should deviate very little from the Hugoniot values. Hence, this 
equation of state should give liighly accurate results. 

Shot 1 was then run on DORF using this refined equation of state, and the results 
again were compared to the experimental values. This time the code gave very satisfying 
results in terms of their deviation from experimental parts. Both DORF calculations are 
compared to both experimental determinations in Fig, 3. Axial pressure plots for selected 
times after impact arc shown in Fig. 4. These pressure plots, as well as all others that 
follow in this report, are derived by assuming that the pressure in a cell represents the 
pressure at the geometric center of that cell. Peak shock wave pressture is taken to be 
the largest cell pressure in a given region of interest, and the shock front position to be 
the geometric center of that cell. 


\ 
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Fig. 4 — Axial preuure ploU for ihot 1 Tor selected times after impact 


It is clear that shock propagation calculations ore extremely sensitive to the equation 
of state used and that great care should be exercised in selecting an equation of state for 
such calculations. An unquestioned use of the Tillotaon formulation would probably have 
resulted in the conclusion that Eulerian codes are not particularly suitt?d for the calculation 
of shock wave propagation problems. It should be pointed out, however, that the Tillotson 
aluminum equation of state is meant to be used over a much wider riuige of material con¬ 
ditions than was necessary for this calculation. Also, its unsatisfactory performance for 
this particular calculation does not necessarily invalidate its use in other calculations. The 
Tillotson equation of state is probably quiU> satisfactory for hypervelocity impact calcu¬ 
lations where the primary interest lies in the computation of projectile or target deformation. 
It applies in these coses because most of the physical processes involved occur at pressures 
much lower than the Hugoniot pressure, wliere the Tillotson form deviates very little from 
experimentally observed values. 

Shot 2 

A 1.27-cm-diameter steel sphere was impacted normally into a 2.B-cm layer of lithium 
hydride, backed by a 4.6-cm layer of Composition B high explosive, at 5.6 km/sec. The 
primary interest was in the peak pressure and pulse shape of the shock wave just after It 
(aib-red the Composition B. Experimentally there were two sources of information: 
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l.27-em>DIAM SUSCL 
SPHERE WirH( 

8.B Km/(*c 


8TEC 

o- 



Pig. & — Configuration for >hot 2 


Table 1 

Pressure Results for Shot 2 


Method of 

Pressure Entering 

Pressure Measurement 

Composition B (kilobars) 

Shock velocity in 
Composition B 

Run distance in 

70 

67 

Composition B 

DORF9 calculation 

70 


(a) measurement of the shock velocity just inside the Composition B yielded a peak 
pressure of 70 kilobars and (b) measurement of the run distance to detonation in the 
Composition B yielded a peak pressure of 67 kilobars, assuming a pulse width of 0.5 usee 
or more. The DOIIFO calculation yielded a peak pressure of 70 kilobars anti a pulse width 
of about. ^ /jsec, (In the I)ORF9 calculation, the Composition B was treated as an inert 
substance. No attempt was made to numerically simulate the detonation process.) The 
shot configuration is shown in Fig. 6, and results are given in Table 1. Shock pressure 
profile.s on axis are shown for various times after impact in Figs. 6 through 9. Again, 
the agreement between calculated and experimental results is quite satisfactory. 

Additional Calculations 


In addition to the above, two other similar calculations were performed without 
experimental counterparts and are referred to as shot 2A and shot 3. The shot configu¬ 
rations are shown in Fig. 10, As can be easily seen, shot 2A has the same configuration 
as shot 2 but with the Composition B replaced by additional lithium hydride. Shot 3 is 
the same as shot 2A, but with a 0.16-cm steel plate in front of the lithium hydride. The 
axial peak pressure decay curves for both shots are shown in Fig. 11. As can be seen, 
the effect of the frontal steel plate in shot 3 becomes damped out at large distances from 
impact. Pressure versus axial distance from impact plots for shot 2A are shown for four 
selected times after impact in Fig. 12. In addition, relative densities for shot 2A at four 
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Fig. 6 ~ Prcuiurc nn uxU v» dUtance iiitu Urg«!t 
for ihut 2 at tlmr 0.88 ittec after impact 
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DISTANCE PROM ORIOINAt PRONT SURFACE OP TAROET (em) 

Fig. B — Pr«Mur« on axit vt dliUnc* Into targvt 
for (hot 2 at time 2.78 iiuc After Impact 



DISTANCE FROM ORtOINAL PRONT SURFACE OF TAROET (cm) 

Pig. 9 — Preiiure on axle vi dUtance into target 
for (hot 2 at time 3.74 vaec after Impact 
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selected times are presented in Fig. 13, The method of display is to plot within a cell a 
number of randomly scattered particles proportional to the relative density in that cell. 
These particles are not carried as part of the calculation, as are the massless “tracer 
particles,” but are created as a means of displaying data only. 



SHOT 3 


Pig. 11 — Peak shock wave pressure 
attenuation curves for shots 2A and 3 


Pig. 10 — Impact configurations for 
shots 2A and 3 
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-1 0 2 3 1 2 3 A 

DISTANCE FROM ORIGINAL FRONT SURFACE OF TARGET DISTANCE FROM ORIGINAL FRONT SURFACE OF TARGET 
(cm) (cm) 



Klg. 12 — Preuure vi distance along axis for shot 2A 









(c) Time > 6.61 Miec 


(d) Time ■ 7.46 liiec 


Filj. 13 — Weighted density plot* for shot 2A at varioui time* after impact 
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Single>Plate Hod Impact Calculationa (DORF) 

The impact of end-oriented metal rods into relatively thin metal plates has been 
extensively studied, with particular interest attached to the length and velocity loss of the 
rods after penetration of the plate. Sets of analytical equations for these quantities have 
been developed which agree quite closely with experimental results. Using the DORF 
code, seven rod impact calculations were made using single-plate targets to determine the 
accuracy with which the code could predict rod length loss and residual velocity. The 
basic configuration and zoning of the problems and tracer particle plots of the impact 
process are presented in Fip. 14 and 16. The equation-of-state constants used in the 
calculations are shown in Table 2 and the numerical results in Table 3. The results are 
quite good, with rod-length-loss calculations differing by no more than about 6% from 
the experimental values. The residual velocity comparison fares even better, but it should 
be kept in mind that both the calculated and experimental velocity losses are extremely 
small. Hence, the residual velocities differ little from the original Impact velocities. 

The DORF calculations made for this study provided the answer to a question that 
has been a subject of discussion for quite some time. The analytical approach used in 
the development of the equations for rod length loss has as one of its hypotheses that 
the extra length loss over and above that which would be expected from steady-state 
incompressible flow theory was caused by the initial shock and rarefaction of impact 
shattering the ftontal portion of the rod. However, the experiments of Christman and 
Oehring (8), where rods were actually x-rayed as they penetrated thick targets, indicate 
that this is not the case. Also, other projectile configuration experiments at NRL indicate 
that the extra rod length loss is not a shock-rarefaction phenomenon. To resolve this 
question, the results of the DORF calculation for an aluminum rod Into an aluminum 
plate wi^ a rod diameter of 1 cm, a plate thickness of 4 cm, and an impact velocity of 
4.6 km/sec are plotted in Fig. 16. The solid line shows the rod length loss as a function 
of rod front position that would be predicted by steady-state incompressible flow theory. 
The DORF data points do not deviate slgnificaotly from this line until the rod front is 
past the original back of the plate. Indicating that the extra rod loss is not a shock- 
rarefaction phenomenon, but is associated with the breakout of the rod from the back 
of the target. 


Low-Velocity Multimaterial Impact (DORF9) 

Two important features of the DORP9 code for Impact calculation purposes are its 
ability to treat up to nine different materials in a problem and Its ability to treat material 
stren^h (using a rigid plastic model). Both of these features were tested by calculating 
the low-velocity (3.048 km/sec) Impact of a steel cylinder into an ablative layer backed 
up with aluminum. The impact configuration is shown in Fig. 17. The equation-of-state 
constants used in the calculation are listed in Table 2. The results are highly strength 
dependent, since the projectile will be deformed and decelerated by the Impact, but will 
do very little actual flowing. Figure 18 compares the calculated projectile configuration 
after impact, using a tracer-particle plot and an actual photograph of the projectile after 
impact. The similarity is very close. Calculated residual length Is 70% versus 68% experi¬ 
mental. Calculated residual velocity is 2.64 km/sec; the experimentally n»asured value is 
2.60 km/sec. 
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ROD 



Pig. 14 — Bulc oonfigurttlon and lotting of kingl«'plit« rod Impact 
oaloulatloni. The thicker and thinner target oaloulatloni were made 
by filling In or deleting material in the void behind the plate, 
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(Page 18 Blank) 



(•) “ 19 waeo (f) Time - 24 /isec 

Pig. 16 - Tracer particle ploU at leleoted timei after Impact 
for a 1.0-cm-diameter aluminum rod 10 cm long impacting a 
1-cm-thick aluminum plate at 4.8 km/ieo 
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7S ORAM STEEL 
CYLINDER -2.9 cm 
LONO « I,IB om RADIUS 


3.048 

-► 



Pig. 17 — Impact configuration for 
loW'Valooity impact 


EXPERIMENTAL PHOTO 



RESIDUAL LENGTH «1,57cm 
RESIDUAL VELOCITY* 2,eOkm/Mc 


DORP 9 CALCULATION 


t 




RESIDUAL LENGTH • l.62sm 
RESIDUAL VELOCITY «2,54km/sec 


Fig. 18 Ilealdual projectile after low-volocity Impact 
D0RF9 code calculutione va experiment 
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CONCLUSIONS 

The experiment-code comparisons that have been presented here obviously reflect 
very favorably on the ability of the codes to predict reality. It should be pointed out, 
however, that only certain specialized phenomena were considered in this study; shock- 
wave pressure attenuation, projectile length loss, deformation and velocity loss, and total 
target penetration. Lateral damage to the target was not considered (the problems were 
zoned finely only near the projectile trajectory axis, leaving coarse zoning in those regions 
off the axis where one might expect the lateral damage to be defined). However, good 
correlation for lateral damage to cadmium plates impacted by cadmium spheres, as well 
as for pressure distribution on a second plate, has previously been reported by one of the 
authors (9). It is felt that the codes are extremely useful tools for hypervelocity impact 
problems, as long as care and an understanding of their weakness are exercised, It is 
highly probable that other fields of endeavor could benefit from their utilization. 
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PH3 


HERE THE DEVIATORIC STRESSES (RIGID PERFECTLY PLASTIC) ARE 
COMPUTED AND THEIR CONTRIBUTIONS TO THE VELOCITIES AND 
ENERGY ARE ACCOUNTED FOR. THESE DEVIATOR STRESSES r; ARE 
FUNCTIONS ONLY OF VELOCITY GRADIENTS, WHERE THE VELOC¬ 
ITIES TO BE USED ARE THOSE FROM PHASE 1. IN TENSOR NOTATION 
THE EQUATIONS, AGAIN DROPPING TRANSPORT TERMS, ARE: 

CONSERVATION OP MOMENTUM: 


au,- 


CONSERVATION OF ENERGY: 





HERE E IS THE TOTAL ENERGY PER GRAM, INTERNAL PLUS KINETIC. 


PH2 


FINALLY, THE TRANSPORT TERMS, THAT WERE TEMPORARILY 
OMITTED IN PHASE 1 AND PHASE 3 ARE SOLVED FOR, AND THE 
MASSES, VELOCITIES, AND SPECIFIC INTERNAL ENERGIES ARE 
INTEGRATED TO TIME t + 

CONSERVATION OF MASS; 


f r • pu - 0 


CONSERVATION OF MOMENTUM: 


dpu 

dt 


-V ■ l^uu) 


CONSERVATION OF ENERGY: 


dp£ 
dt " 


V • fjxiE 


t 

RECYCLE 
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